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In the third quarterly report, a brief progress report on the M.S.U. laboratory 
facilities will be presented in addition to the use of dry processin,T and alternate 
dielectrics. Emphasis will be placed on the two-dimensional MOSFET simulator. 

The two dimensional modeling program has been written for the simulation of 
short channel MOSFETs with nonuniform substrate doping. A key simplifying assumption I 
used is that the majority carriers can be represented by a sheet charge at the silicon , 
dioxide-silicon interface. In solving current continuity equation, the program does j 
not converge. However, solving the 2-dimensional Poisson equation for the potential | 
distribution has been achieved. The status of other 2-D MOSFET simulation programs | 
are summarized . ! 
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1. MSU UBdl^KSr fACILIfXBS 

negotiations are in pregrass vith two saaieoMuctor coapaaias 
to donate surplus production aquipaent which Include - (1) epi 
reactor and generator, (2) ion punp Vacuua systaa, (3) vapox reactors, 
and (4) aiscellaneous equipaent for the hybrid laboratory. 

- X 

II. DRY PROCESSING AND ALTERNATE DIELECTRICS 

- The exploratory processing involving dry cheaisrry required 
the use of plasna etching and plasaa deposition equipsient not 
available at either the Marshall Space Flight Center or Mississippi 
State University. Consequently, the logistics involved in pro- 
cessing the wafers have been difficult. 

Four types of insulators were to be investigated: 

1) 5Z phosphorus doped silox, 

2) plasma nitride, 

3) the GAF version of polyimlde, and 

4) the Hitachi PIQ. 

Whenever possible, dry chemistry processing was to be compared with 
"standard" wet chemistry. The wafers were processed through the 
first metal mask at Marshall Space Flight Center in standard 
fashion. The polyimlde material was processed at American Micro- 
systems, Inc. (AMI) using the standard recommend procedures by 


the MRttf Actu feea* The eltrlil eilox iMM ii^ilted et tM 



plewi rMtttor. TIim iiet eheaifteiT pet^omed at 

AKt ttstot ^'a ataadard ptoeaaaiag. The plaaaa etching vaa 
perfotned fay Davla and Wilder on one of their production parallel 
plate plaana ayaten using their proprietary processes. 

The vet chentstry split involving sllox vas lost when the 
standard "pad” etch attacked the aetal coapletely. The aaae 
batch etcttant continued to provide good results on Mfl aaterial. 
Since the first layer a«tal at NASA was deposited via sputtering 
and with many different impurities » conpared with approximately 
1$ Si:Al deposited with electron beam, the difference in the 
metal film is believed to be the source of the problem. 

The plasma nitride material and silox naterila, deposited by 
AMI and plasma etched by Davis & Wilder, vas processed without 
any problems other than long delays. This material plus the 
polylmide material has finally been received for the final 
processing here. Although, no electrical results are available 
at the present, the microscopic examination displays the expected 
sharper edges of the dry processed material vs. the standard wet 
chemistry. 

In the meantime, the results of the in situ back sputtering 
experiment were obtained. As vas expected, the in situ back 
sputtering split gave much Improved yields over that for the 
standard material - the detail results being published elsewhere. 

The results of these experiments points to the requirement to 
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H«t7 «ff«6tiv«ly r«M9v« th« AljP^ prior m tbo xifo i ia MtaUsttlSR 

or i dnia lre I f for—t lwi durli^jmfvr pro£Mfiai-Je-«i£tiiit.rlio 

• ^ 

order of 30 A idiich a good sititoting traataent egtl tadoea. Tha 
In aitu elaaning forgivaa prior procaaaing iditeh tanda to g;ro«r 
tha Al^Oj such aa idMtoraslst raMoal and buf farad IF atehli^ 
of oxide. 

tte dry ptoceaslng haa the potential in tlw fti^re of 




reaolving aaaller via sizes and* sisMltaneously) avoids the 
haraful B2O to accelerate the altoiinaBi oxide (and hydrate) 
fomation. Although dry processing to reaove the photoresist is 
routinely done one shotild expect Al20^ fomation due to the 
presence of active oxygen in the "ashing" process. Some parallel 
plate reactors can be modified to Introduce Ar at lov pressure 
and effectively as in a back sputtering mode. This latter mode 
could then be used as a clean up step prior to second metal 
deposition - although not in situ. Undoubtedly, dry processing 
is going to be extremely important in future double layer metal 
development work. 
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III. TWn-DIMENSIONAT. MOSFRT STMUUTTilN PROGRAM 


A. Long Tent Objective 

The basic long term objective is to develop a two-dimensional 
computer program which derives the channel current in steady-state 
from given values of potentialr source, drain, gate, and substrate. 
Also given are physical parameters such as gate length, dielectric 
thickness, doping levels, etc. The substrate doping level should 
be variable in two dimensions and the substrate potential is 
independent of the source potential in order to include body-bias 
effects:. The program is intended to model MOSFETs with short channels 
with particular emphasis on subthreshold and punch-through conduction. 

Secondary objective is to utilize techniques in the computer 
program to minimize computer run time while still maintaining reasona- 
ble simulation accuracy. 


B. Introduction 

At the beginning of this project no programs were available 
which satisfied the specified objectives. There has been consid- 
erable work in two-dimensional modelling of JFETs including the 
works of Kennedy and O'Brien [ll and Martin Reiser [2]. The 


iffiSFBT Hii itw ilittlat ioaV 

•.g.» Barnaa t3)* Varioua ao4aa of ttOSFBT o^raci^ hava baan 
alaulated or nodeled using twoHIlaatMlonal teahaitittas, a.g. 
Schroader and Muller [4] (saturation), Barron (51 (low level 
currants) , Araatrong e al (6) (pinch-off), Kennady 17] (affe«*t8 
of ionizing radiation), Kennedy and MSrley [8] (saturaticn) , 

Mock [9], Motta [10], El-Hanay and Boothroyd [11]. Although 
these works all concribute to the knowledge required in reaching 
the objective of this task, only the effort by Mock is sufficiently 
general to meet the objectives — and it is unavailable for public 
distribution. 

The work by Barnes [3] described the application of the 
finite-element method to the analysis of the MESFET. 

Further, the finite-element method has been used for some time 
in solving problems in mechanics and elasticity; however, it has 
only recently been applied to semiconduction problems. This method 
has the power to treat some problems, such as eigen-value problems, 
for which the finite-difference swthod is awkward if at all 
applicable. It can also be applied to the solution of field 
distributions governed by partial differential equations, and one 
of the most attractive features as compared to the finite-element 
method is purported to be the ease of treating non-rectangular 
geometries and irregular boundaries. For example, the geometry 
of the VMOS structure could be accomodated. It was decided to 
further investigate this technique. 
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Tlii Jitillt of thfi irralfiiftira iro in ttM 

In nu—iry, Clio «vnltintion r«Mlto4 In mm nteptlciM ebnt CIm 
f ittitt olement M€hod tmuld b« nuf f Iclontly of fcotlvo for aMi- 
conductor probloan to Justify tho effort. A M>re recent pepe<‘ hn* 
reinforced this attitude. This paper indicates thet the proper 
formulations of the semiconductor problem for the finite element 
approach remain to be demonstrated, and, in agreement with out 
observations, point out that the application of Galerkin's method 
is subject to skepticism. Therefore, it was concluded that the 
further work should be based upon the finite-difference method 
which we have used before although the finite-element method is 
intriguing and may be further developed in the future. 

Key in the development of this task is the proper selection 
of assimiptions. Certain assumptions are required to focus the 
effort and facilitate convergence expeditiously and, therefore, 
minimize computer time. On the other hand, too many assumptions 
lead to loss in utility. 

The basic assumptions follow: 

1) Recombination/generation is negligible. 

2) Hole conduction is negligible. 

3) The conventional expression for Poisson's equation is 
valid [see next section]. 

4) The mobile malority carriers are included in an 
infinitesimally thick layer of charge at the S 1 -S 102 interface. 

5) The current flow in the channel is then described by a 
one-dimensional equation. 
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S ■ density of majority earricrs in the ehaimclt 

D • majority carrier diffoaiem constant • 

U • channel carrier mobility, and 

• tangential electric field component. 

The first three sections are common in the stated literature 
and restrict the model from covering phe n om en a associated vith 
breakdown and bulk generated leadage currents. The last assumption 
is key and separates this work seriously from that of Ibck although 
Kennedy (7] used a similar approach for the sake of simplification. 

Key also arc the normal assuaqitions which are not made and 
tHilch are conaionly utilized: 

1) the gradual channel approximation, 

2) constant substrate doping, and 

3) zero substrate bias. 

Although a constant effective channel mobility is utilized 
during the initial development of the program, it is not Intended 
to remain a basic assumption. 

Assuming step junctions at the source (or drain), the potential 
in the depletion region and c/'te is defined by Poisson's equation 
with apprcpriate boundary conditions (Figure 1): 

1) fixed potential along the source (or drain) boundary 
(BCD and UK), 

2) fixed potential deep within the bulk substrate (AU), 
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3) Mvo horitontal eeapoiwat of oloetrlc field wider the 
•owree end drain and far fron the ehannal tAS and Kt>« 

4) aero horlaontal conponant of el.*c ^Ic field in the oxide 
far from the channel (EF and m) and 

5) constant potential on the gate electrode (FG). 

For convenience the substrate in the basic progran is defined 
as having zero potential. The source, drain, and gate potentials 
are heir respective applied potentials plus their irark functions 


with respect 

to the substrate 9 

i.e, 

' ♦ 


¥ « 


+ 9 


source 

S 

Bi * 



V 

+ V . 


drain » 

D 

Bl ’ 


^ * 


■f ^ ♦ 


gate 

G 

ms 

For the SMt 

example the potential 

is always non^negatlve* 


The gate dielectric is assuned to be ideal with infinite 
resistivity. The charge density in the dielectric is assuowd to 
be zero except for a sheet charge representing the silicon-silicon 
dioxide interface charge Q which is arbitrarily assumed to be 

o 

50 A within the oxide at one grid spacing above the channel grid. 



DISTANCES 
AB = YDIST - YDTF 


uO-ORDINATES 
A = (],]) 


CD = JI - YDIF 


B = 0 , JKl) 


EF = HG = TOX 


C = (IK.JKI) 


= JK = XDTF 


D = (IK.JMAX) 


AL = XDIST 


F = (IK5,JL) 


OTHER RELATIONSHIPS 


IKS » lK-5, IL5 = IL +5 
JL^* - JL-1, JK = JKl-1 


JMAXl = JMAX +1 


JMAX<!» = JMAX-1 


G = (IL5.JL) 

I = (TL.JMAX) 
K = (IMAX,JK1) 
L = (IMAX.l) 


FiRure 1. MOSKET Cross- Sect ion 
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Key to the slnulatlcm of short channel devices Is solving the 
tvo-dl9K!n8lonai Poisson equation vittK>ut asking t^ecessary 
approximations* The finite^difference formulation for Poisson’s 
equation is given in the next section. The majority carrier charge 
is neglected in this formulation except for the sheet charge 
representation at the oxide interface. The intent is to ignore the 
charge where it is negligible cosq>ared with the ionized doping 
charge — which is a good assumption everywhere except in the channel. 

One can then solve Poisson’s equation for the two-dimensional 
potential distribution more easily and faster than if the coupled 
majority carrier current equation is included and the majority 
carrier density variable is included for every mesh node. 

After finding the potential distribution one can» after the 
fact, explicitly solve for the majority carrier densities and 
current in order to predict punch-thru and subthreshold conduction 
behavior. This approach is aimed at minimizing computation time 
whi?e retaining the essential relationships to maintain sufficient 
accuracy for modeling purposes* One should expect a compromize in 
accuracy, under certain conditions, when compared with the more 
general approach of solving the two-dimensional conduction equation. 

C. The Poisson Equation, Difference Form 
The Poisson equation in different form is 

» - r/t (1) 
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f ■ _ electrostat^^ potential^ ^ 

p “ charxe density, and 

e ■ peraittlvlty of the Material. 

For acceptor doped material tite charsa densl^ Is 

“ q(p»idl^) , depletion, (2) 

where 

q ■ charge density, 

p “ hole density, 

n = electron density, and 

= ionized acceptor density. 

The hole density is defined by 

p « exp(-u) (3) 

where 

u = qV/kT , the noraalized 
potential, and 

u = 0 in the charge neutral bulk. 

Consequently, in the bulk substrate and in the depletion layer where 
the electron density can be neglected, 

P = - q Il-cxp(-u)] . ( 4 ) 

In the channel region the assumption is made that the electrons 
exists only In a sheet form of infinitestimal thickness at the 
silicon dioxide-silicon interface. Since the hole density is 

Insignificant at inversion, the charge density in the channel is 

Pg --qCNaS+ns) t channel, (5) 
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where 


R^g- is- the effective sheet cherge due to- the I o nise d" 

scceptors within the grid spacing defined by the interface 
grid, and n^ is the sheet charge of electrons in the 
channel. 

Near the oxide-silicon interface there are positive interface 
charge ~ of density N^g , consequently. 


qNgg , interface. 


( 6 ) 


with the location being assueed to be 50 A inside of the oxide. 
Otherwise in the oxide the charge is assumed neutral. 


0 , oxide. 


(7) 


To solve the Poisson equation numerically, the left hand side 
can be expanded in difference form as 


•> 2 
V *!■ « V„ V U 

T 


( 8 ) 


■ "t “ "i-i.j * <= “i,3 




where 

A 

B 

D 


kT/q , 

h^)J 


-1 


2[Wi_i(Wi_i+ , 


2[Wi(Wi_i + Wj)l 


-1 
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B - h4>r . 




-(A + B D + B) 


vith tlie grid systen defined in Figure 2. 

Ue^og the charge equation (14) the coaqjlete Poisaon equation 
in finite-difference fom is written for node (i»i) as 




+ E(i,j)U, 




a 


[ 1-exp^-U. .)] 


eV» “ ™ “i.j- 


(9) 


i»j^l ‘'’T 

The above nonlinear equation can be linearized by defining a function 
F such that 


^^“i,j-l’ Vl.j* ®i,J» ”i+l,j* “i,j+l^ 


i.j 


^(i.j) '^i,j-i ®(i,j) + C(i.j) Uj j 


^ ®(ij) “i+l,j * ”i,j+l 


- (qN^/eV^) (1-exp - 0 


( 10 ) 


Defining a set of reference value of U's and F such that 


uy, 4 » «tc. and fJ , one can expand F around the 

i»j“l l”i*j l»j i»J 


reference potentials: 
















‘•I.1-1 • Cj-I - «I.J-1 • ‘“' 


The above exjureselone def iitee e syeteB of equatloM In AU as 


^.j ■ ‘"i.j 



where 


A - A, B = B. D - D, E « E, and 


C - C - (qNa/kT) exp (-uj j) 

with A through E being defined inequation ( B )• The function F is 
theoretically zero; therefore, by setting ^ * 0 , a system of 
linear equations in AU is obtained. For the node l,j 




* "'"'m.j ^ “"i.JH-i ■ - ■'i.j 

where F^ Is defined in equation (10 ). 

^ »J 

By applying the method of lines, the linearized equations are 
resolved into a tridiagonal matrix form as 




-"i.d -'*"1-1,1 -“'"'i+i.j 


where AU are taken as the roost recent updated values of AU 


each iteration, but still keeping the U’s fixed until convergence 
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it aehiftvtd. above aquation (U ) ratal j^ tt^ old 
9 ^ AO along tha x»axla and calculataa tha nev v alua a of Att_ 
alof^ tha (^vlously , aquation (14 ) can ba writtM for 

tha tnvarae caaa. By altamatlng batwaan tha vertical and 
horliotttal foma of the tridiagonal matrix, tha convarganca can 
be axpaditad. 


Special conaldaratloh muat ba givan to Poisson's equation at 
the silicon-silicon dioxide interfaca. Figure ( 3 ) Illustrates the 
channel sheet charge S in charge per unit area, the acceptor 
charge N in charge per unit volume, and the discontinuity in 

o 

permittivity. Applying Gauss's law in one dimension with two 
surfaces at the midway points between C & C 1 and C - 1 & C 
leads to the expression 


^0x^1/ 2 


^Sv-1/2 


positive charge enclosed 


(15) 


M h , 

- q(S ) 


Multiplying through by 




leads to 


h +h - 

C C“1 


U . U , U 
ox ( c»l- c) c-1- c 

?i c c-1 




S+- 


N h ,) 
a c-1' 






(16) 


which is the one dimension form of Poisson's Law 

dy^ Si 


( 17 ) 
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2«|(8 4 * 

‘^elf " " V~ (h 

If C" ©« t 


Therefore, the finite difference form of the eqeation is dhtaitied 
by modifying the folloirifig terma te eqnetlon ( 9 >| 




C - -(A + B + D + B) , and 


!eff . ^ (S ^ N^^.,/2) 

®Si ^T*^Si ^**c'*‘ 

[Note: y J for the channel node in the program. ] 

iBax 

Similarly the linearized version is modified by 


« E^, Cj - Cj, and 


F" = A.U. , + B.U, , + C.U. 

i,c i i,c-l i i-l,c 1 l,c 


* ^ ".ff '“^Sl 


where the subscript i has been Inserted to denote the coefficients 
for mesh node (i»c)e 

Tn the region of Si - SIO2 interface where the charges 
are assumed to exists » one mesh node is assumed inside the oxide, 
the coefficients are 


A * A, B * B, C * C, D » D, E * E, and 
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' - ^I.c^ «*l-l.cHa^ 


2Q-S . 

+ — f>» + h^^fl 
e V* € e+1 

ox T 


SlBllarly , in the oxide where p ■ 0 , the coefficieata 


A ■ A, B ■ B, C ■ C, 0 • D» E - E, end 
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D. The Channel Conduction Equation* 

Difference Fora 

As noted previously* a significant assunption utilized to 
facilitate convergence, is that the channel charge is modeled as 
a sheet charge at the silicon-silicon dioxide interface. The 
channel current is assumed to obey the single dimensional form of 
the conventional three dimensional equation* i.e.* 


IS 


( 22 ) 


where 


S « density of electrons in channel (Iper unit area). 


D = electron diffusion constant, 

n 

■ electron mobility, normally field dependent, and 
E ■ tangential electric field component. 


The above equation can be written in normalized form utilizing 


E^ ■ — and (Einstein's relationship): 


19 - 


Thffl 




I - (s-vsiu^ 

s n 


(2S) 

(24) 


wlicrc - J^(qV^)“^ . s'- ^ H 


as 


»u 


AC «iy point X In the chwuMl th« t<tAor«il oolutioa of S la 

8(x) - axp / I L ( IT ®*P /* ® *t + cj (25) 

* 1-1 ( * 1-1 “n * 1-1 ) 

which reduces tc Che specific solution for S(x), given 
S at X ■ X. 


‘ 1-1 * 

-/ N -U e X **“ 

S(x)e - ^ 


-u(0 

7(0 




i-1 


Rewriting the above leads to in terms of ^ 


( 26 ) 


S^exp(-U^)- Sj^_j^exp(-Uj_j) 




ex£ 


ziLCU 


AS 


*i-l 


<J V- 


( 27 ) 


Similarly, integrating between and leads to 


S^^lC-xp.(l)i^,)- e,p-(U^)_ 


i+1 


exp-U(f ) d ' 
m(0 


qv, 


( 28 ) 


Assuming current continuity, i.e., J ■■ constant, leads to 

n 








I 


’i-1 


(») 


+ («ip - u. 


‘l+l 


«^p(U(Ci dC 

wCO I 


‘ 1-1 


f* 


- (exp-U 


141 


1 exgtulsDd CJ s . 0 

1.(0 ( ®14i 


‘i-1 


The above equation is of the form 


*i®l-l **■ ’'"i^i ■‘‘ ®i®141 ■ ® 


and define a tridiagonal linear system of equations. 

Equation 27 should reduce to a fwnlllar fora as ‘ 

Let S. ,» S , S.»S 4 .\s, U. ,■ U , U.» U 4AU . Then 27 becomes 
i-1 o 1 o 1-1 o 1 o 


-I’ -U 

(S 4 AS)c ®(1-AU)-S e ° J 
o o _ n 


-1 -U 

u e o Ax 


qv. 


(30) 


which reduces to 


„ jAfS - AUl 
I Ax " ®o Ax f 


(31) 


which displays the familiar diffusion and drift terms in the limit. 
The integral form in ( 27 ) has the potential of being more accurate 
with larger grid spacing. 

Another observation one may note is that the coefficients 
^i * ^1 ' ^i scaled, or normalized, to values less dependent 
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^ potMtUl by UMi W «^<^|) 

(1<WM ~~ 



•j • - mtpiV^ 




U(0«syft!4-V<c)]d € 


.1 

**1 " ** *Cr Mi 

*t-l 



*^1 






* i"^(C)eM IUj-U(€)J ie . (32) 
* 1-1 


And, of course, if the aobility le ast«aed conetmt, then it can be 

renoved by multiplying through by u . 

Noting chat the integrand varies exponentially with potential 

whereas the potential is expected to vary with x no worse than 

the third power, one can approximate the integral by assuming 


- r ■ a(x - X ) (33) 

o o 

dn 

and evaluate the slope ^ ^ ^be lowest value of potential 

over the range of Integration, In this manner Che neglected higher 
order terms contribute very little error since the integrand falls 
off exponentially with potential. Changing variables for Che 
integration of , as an example, leads to 


f'**' £iia 


exp Cl’ j * I’l dl? 


(34) 


1-1 


a 


dU 

dx 


Further, 


where 


(35) 


{exp(Uj- - eKp(P^-tlj^j^)) 

vith a and u evaluated at the lowest value of tJ which for a > tX 
la . Sisllarly, 

“ (aul'^lexpdfj-r^^j^) - 11 exp(Uj- 

-3 

Ci = ■ «*P^V ”i-l^* exp(U^- 

“1 

bi = [aplj^_j^texp(Uj- \ 1 ^_^) - erp(0^- (36) 

where the subscript on (ap] denotes the node for evaluating a 
and p. assuming a > 0- . 

The boundary conditions for this systea of equations 
association with electron conduction in the channel must be defined 
in terms of the source and drain corner mesh nodes for the MOSFET 
operation in the active region. Figure 4 illustrates the applica- 
tion of Gauss's Law at the source corner of the channel. Mesh 
nodes along the boundary of the source are common poential, therefore, 
the tangential electric field terms along the boundary of the source 
are zero. Since the electric field along the channel F.^ for 
operations of interest - at least at this source corner, the surface 
charge is modelled in a single dimensional manner, i.e. , 

e,,E«--qS (37) 

or, in terms of potential. 


V-e (F 
*1* ' «■ 




* qS 


(38) 
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All alternate scheae involves the solution of the one 
dlaensional Poisson equation: 

- ^ I « V ♦«> - V ''bi> I - * "b («> 

OX » ) 

ubere V • the applied bias to the gate 
G 

4 ” the gate work function relative to the substrate* 

QS 

V * the applied channel potential , 

v» 

V,. , « th 2 work function of the drain/source, 

Bi 

t * oxide thickness y 
ox 

= oxide permitivity, 

Q = silicon dioxide-silicon interface charge, 

q = electronic charge, and 
Q = the bulk charge density (per unit area) 

» 0 for drain and source. 

The expression can be defined in terms of the vertical electric 

field in the channel region: 

Q = eV (^C"l " ^g) - q Na }lSzL 
c-1 

where the subscript c refers to the channel node. 

E. Channel Coupled Equations 

Although the channel conduction equation f29) is linear in S, 
it is nonlinear in potential. The simultaneous solution of the 
conduction equation and the Poisson equation along the channel 
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requires linearising ( 29). 




*1-1 + ‘‘i ®i ♦ 'i *1+1 '*•” V *1’ 't 


are defined in (36 ). Since the teraa fan] vary acre slouly 
than the exponential texas* the tens is asstaed constant for ai^ 
particular iteration, using the last knoim values of u'a for Its 
evaluation. Then, the linearized conduction equaticm becoaes 


1 1 

*ci ■ ■ - * 1 * 1 - 1 -<">I- 1*1 “>*‘V"i-i> 


+ CjS^^jexp(U^-Ui i) U-exp(Uj-Oj_^)] 


-1 


B * * a* 
cl 1 


r-F*” 


.- 1 . 


3Ui 


- 1 . 


‘"i^l+l^^ “ exp(U^-U^_^)] ^11-2 exp(U^-U^_j)J 


"cl “ ‘’l ‘ 


r. 

E . “ £1 

Cl 




-1 


+ bjSj exp(Uj-U^^^)Iexp(U^-U^_j^) - exp(U^-U^^j) J 


-1 


- ‘■i*i+i ’ 


n . * c. , and 
ci 1 


= F*’. + A .AU. , + B .AS. , + C ^AU, + D .AS. (40) 

Cl cl Cl 1-1 cl 1-1 cl i Cl 1 


E .,AU_, + C .AU,^, =0 

cl 1+1 cl 1+1 - 
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The Poisson equation has been previously derived, equations 
QS ) nod (l9 )> with the coefficients being defined below for the 
new system of equations: 


Wl.c ^i^«ic “i^^i ^ ®i^“i+l.c 


where 


= -A.AU. - K.AU. . , - f' = 

i i»c-l i lgC+1 i,c ic 


“l- ♦ V-l” ’ 


■ A.^ from (17 ) , 

from (17 ) , 

from ( 17 ) , 

from ( 17 ) , 

= Cj from ( 17 ) , 


F = as def Ined in ( 28 ) 

l.c 


Simultaneously solving for the channel potential and Che 
change charge requires Che simultaneous solution of the above two 
equations (40 ) and ( 41 ) which results in the following system 
of equation in matrix form: 


C H D 0 0 0 


C D E G 0 0 . 


B 0 C H D 0 


A B C D E G 0 . 


OOBOCHDO. 

OOABCDEGO 

////// 


N-l,c 

^N-l.c 



TiM Matrix is a 6 diagonal Matrix* tHa fiixibar of variablas is ? 
twico-tiia ntad>sr of ^id 

the boundat^r nodes at the source or draitt* TlM solution 
conveniently solved by using a Guassian el Ininatlon technique. 

F*. Overall Flow Chart 

In order to expedite convergence, initially, the gradual 
channel approxiaation is assuned and the potential is calculated 
along the channel. This analytical approach divides the device 
cross-section into two regions, the channel, silicon substrate 
region and the oxide region. 

In each of the two regions an initial approxinatlon of the 
potential distribution is obtained by solving the one-dlnensional 
form of Poisson's equation utilizing the given boundary potentials 
Including the channel potential derived in the previous step. The 
one-dimensional equation is equivalent to setting dU/dx > 0 . 

The potential distribution is found along one column (vertical) 
simultaneously. By sweeping in the x direction, i.e., repeat 
Che solution on another column, the potential distribution for 
the entire region(s) is obtained. 

Next the solution to the two dimensional P. E. is obtained 
in each region. In the oxide the potential along a column is 
obtained with the region being swept in the x-dlrection. In the 
silicon convergence is expedited by using the alternating direction 
implicit (ADI) scheme. First, the potential in the row below the 
channel (JMAXO) iq obtained. Next the potential distributions for 
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each tov hetween the aourcc and dtain ate obtained by aeeeytiig la 
the ne^ti^ y^^ aolution to toe iii la obtelned» 

the solution to the coluan t • 1 is obtained The colunp solutions 
ate obtained by sweeping in the positive x^ditection* This ADT 
schene is repeated \mtil satlsfactoxy convergence is obtained. 

With the t«o**diaien8lonal potential distribution obtained for 
both regions based on the gradual channel assuq>tion, a new 
estimate of the channel potential and the channel charge is obtained 
by solving simultaneously the P. E. and current equation along the 
channel row. The 2-D P. E. is then solved along the rows between 
the source and drain by sweeping y in the negative direction. Next 
the columns solutions are obtained sweeping the x-direction for 
the entire regions — thus crossing the channel in the appropriate 
locations (1 = I K + 1 to I L-1 ) . This ADI scheme starting at the 
simultaneous solution of the P. E. and the current equation along 
the channel row is repeated until satisfactory convergence is 
obtained. 

The channel current is obtained from an explicit solution to 
the channel current equation after the channel charge and potential 
are found. 


G. Availability of 2-D MOSFET Simulators 

The first 2-D simulator program which could provide I-V 
char.'icterist It'S was developed by Mock (9] and Kennedy (&] at 
IBM in the 1970 to 1972 timeframe. The program, based on finite 
difference techniques, is proprietary to IBM. However, a modified 
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•ubstrate bias. 

Rennady at the Applied Electronics Research Corporation has 
■odified the Sutherland program to include provisions for 

1) 1-D substrate ioq>urity profile, 

2) nobility dependence on gate voltage, and 

3} substrate bias 

and is now licensing the progran (NEMOS) for $10,000 for one year 
and $7,500 for each succeeding year. 

In addition, Hitachi has developed a finite element version 
with provisions, at least, for substrates with 1*0 Inpurity profiles 
Dr. Asal of Hitachi has indicated the possibility of licensing the 
program to Universities, govemnent and industry in the United 
States. The Hitachi program called CADDETvas developed in 1977-1978. 

At least: one other 2-D simulation program is under development. 
Jim Meindl at Stanford Is under contract with ARPA to continue the 
work by Sutherland. 

At this writing the autiior has been unable to determine when 
the program under development at Stanford will be available for 
distribution. 
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